(N 

o 

(N 
S3 

a 



c3 



o3 



S3 
O 

o 



> 
in 

o 

(N 



Catalytic Growth of N-doped MgO on Mo(OOl) 
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A simple pathway to grow thin hims of N-doped MgO (MgO:N), which has been found experimentally to be 
a ferromagnetic d° insulator, is presented. It relies on the catalytic properties of a Mo(001) substrate using 
growth of Mg in a mixed atmosphere of O2 and N2 . Scanning tunneling spectroscopy reveals that the films 
are insulating and exhibit an N-induced state slightly below the conduction band minimum. 



Recently, it has been found that MgO:N-films grown 
by molecular beam epitaxy (MBE) exhibit ferromag- 
netism after being annealed at 1020 K. 1 The optimized 
N-concentration was 2.2 % exhibiting coercive fields as 
large as 60 mT at T = 10 K and magnetic moments 
per N-atom of 0.3 hb barely reducing up to room 
temperature. A Curie temperature Tq — 550 K has been 
extrapolated and indications that N is incorporated 
substitutionally on the O-site have been deduced from 
core level spectroscopy. Moreover, independent studies 
of N implantation (80 keV) into MgO led to a hysteresis 
with a coercive field of 30 mT at 300 K. 2 This raises hope 
that reliable d° ferromagnetism avoiding d-metals can 
be realized in MgO:N at T = 300 K. Such magnetism 
without d-orbitals has previously been found in thin 
films of undoped oxides 3 including MgO 4 or defective 
carbon systems 5 , however, with limited control since 
relying on defects. ZnO with sp-type dopants as C, N, 
B, Li, Na, Mg, Al, and Ga shows ferromagnetic signals, 
too, 6 ' 7 but, likely, Zn or O vacancies and Zn d-orbitals 
are involved in the magnetic coupling 6,8 . 
The d° ferromagnetism has first been proposed the- 
oretically relying on the double exchange mechanism 
in narrow impurity bands. 9 But the high Tq proposed 
originally has been challenged by going beyond the 
mean-field approximation 10 or by considering correlation 
effects, 11 " 13 . Partly, even the absence of ferromagnetism 
was found. 11 This renders the high Tq observed exper- 
imentally into obvious disagreement to current theory 
asking for more detailed studies . 

MgO:N films, in addition, exhibit bipolar resistive 
switching behaviour 14 prior to annealing. Resistance 
contrasts as large as 4 orders of magnitude, switching 
currents as low as 100 nA and switching times into both 
states below 10 ns have been obtained. 1 This makes 
MgO:N also interesting for nonvolatile memories. 
However, the incorporation of N into MgO is difficult due 
to the strongly endothermal incorporation of N atoms 
with respect to N2 (energy cost per N atom: 10 eV) 15 . 
It requires, e.g., atomic beams of N and O produced by a 
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FIG. 1. (Color online) a) STM image of 7 ML MgO0.973N0.027 
film on Mo(001) with thicker islands (bright areas) on top of 
a wetting layer of 1-2 ML; (50 x 50) nm 2 , U = 3 V, I = 0.5 nA. 
b) Line profile along the black line marked in (a) with height 
of the island above the wetting layer marked. 



high-frequency ion plasma source 1 or N + implantation 2 . 
Here, we demonstrate a simplified pathway using the cat- 
alytic abilities of a Mo(001) substrate. Thus, we es- 
tablish a model system of MgO:N for surface science. 
Mo(001) is chosen since thin MgO films of high quality 
can be grown epitaxially due to the relatively small lat- 
tice mismatch of 6 %. 16-18 Moreover, catalytic properties 
of Mo with respect to N2 are known, as, e.g., for the 
nitrogenase within bacteria using molybdenum enzymes 
as catalyst. 19 Catalytic N2 dissociation on surfaces has 
been induced successfully for the electronically similar 
W(001), where growth properties of MgO are, however, 
unknown. 20-22 We have grown thin films of MgO:N on 
Mo(001) with thicknesses up to ten monolayer (ML) at 
optimal doping. Scanning tunneling spectroscopy (STS) 
revealed that the Fermi level is well within the band gap 
indicating insulating behaviour of the film. Moreover, an 
unoccupied state close to the conduction band is found, 
which is not present in pure MgO films. 
The experiments are performed in ultra-high vacuum at 
a base pressure of p = 5T0 -11 mbar. Firstly, the Mo(100) 
crystal was cleaned by cyclically annealing within O2 at 
initial pressure po 2 — 5 • 10- 7 mbar and 1400 K, followed 
by flashing to 2300 K. 23 After every cycle, po 2 is slightly 
reduced. The MgO:N films are prepared by molecular 
beam epitaxy of magnesium at po 2 = 1 ■ 10~ 7 mbar and 
N2 pressure p-^ 2 = 5 • 10 -6 mbar. The deposition tem- 
perature Td is 300 K, if not given explicitly. The depo- 
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FIG. 2. (Color online) a) Differential AES spectra of N and 
O for x — 3.1% with the corresponding peak to peak height 
Yi marked, b) Nitrogen concentration £ of 7 monolayer thick 
MgO:N films as a function of deposition temperature Td. 



sition rate of Mg controlled by a quartz micro-balance 
was 0.5 ML/min. After MgO:N deposition, the samples 
are annealed at HOOK for lOmin. Figure 1(a) shows 
a scanning tunneling microscopy (STM) image of 7 ML 
MgO0.973N0.027- It exhibits islands on top of an MgO:N 
wetting layer. From comparison of the coverage de- 
termined by the quartz balance and the volume of the 
MgO:N islands, we estimate the wetting layer thickness 
to 1-2 ML. We checked crystalline quality and chemical 
purity of Mo(100) and the MgO:N films by low energy 
electron diffraction (LEED) and Auger electron spec- 
troscopy (AES). The nitrogen concentration x within the 
MgO films was determined by AES at primary electron 
energy of 1 kcV being sensitive to the upper 2 ran (9 to 
10 ML) only 24 . The atomic concentration Xi of element 
i within homogeneous films is calculated from the AES 
peak height Yi shown in Fig. 2(a) according to: 25 
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(1) 



where Si denotes the normalized sensitivity factor and a 
sums over all relevant elements. 

The N concentration x depends on deposition temper- 
ature T D as shown in Fig.2(b) for a 7 ML MgO:N film. 
Up to x = 6% is achieved at Td = 850 K indicating a 
more effective dissociation of N2 at higher Tq. Notice 
that x = 6 % at 7 ML is still far below the amount of 
N expected from a full N coverage of Mo(001). Next, we 
prepared MgO:N films with a thickness of up to 100 ML 
at T D = 900 K. The nitrogen amount within the MgO:N 
decreased dramatically with film thickness as shown in 
Fig. 3. A 60 ML thick MgO:N film contained only 0.3% 
nitrogen in comparison to 3.2% at a thickness of 7 ML. 
Assuming homogeneous distribution of N, this implies 
that the amount of N in both films is roughly the same 
evidencing that N 2 is dissociated on the Mo(100) surface 
only. We assume that during MgO growth or anneal- 
ing, the atomic N from the surface is incorporated into 
the MgO film. To support this scenario, we fit the data 
assuming a constant areal concentration n of of N atoms: 
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FIG. 3. (Color online) a) Nitrogen concentration iasa func- 
tion of monolayers d; different symbols mark different prepa- 
ration methods; fit curve assumes a constant amount of nitro- 
gen atoms (see text), (b) Sketch of N2 incorporation: left: N2 
dissociation on Mo(001); right: incorporation of N into MgO 
leaving the amount of N independent of MgO thickness. 



where d denotes the number of MgO ML and n repre- 
sents the concentration of N atoms with respect to the 
sum of N and O atoms. If we assume that all N atoms 
are substitutionally incorporated into the first monolayer 
of MgO:N, the fit curve plotted in figure 3(a) shows ex- 
cellent agreement with the measured data points using 
11 = 21.6%. Assuming interstitial impurities, i.e. dumb- 
bells of NO 15 , equation (2) has to be slightly modified 
and results in n = 18.2 %. Nearly the same n can be 
achieved by another preparation method: Prior to the 
deposition of Mg in a pure O2 environment, we exposed 
the Mo crystal to N 2 (p = 5 • 10" 6 mbar) at 300 K for 10 
min. Afterwards, we grow MgO without N2 at T = 900 K 
leading to very similar N concentrations as shown in Fig. 
3(a). Thus, obviously, the dissociation of N2 takes place 
at the Mo(001) only. Finally, a third preparation has 
been performed: 10 ML of pristine MgO were grown first 
at p Na < 10" 11 mbar. Subsequently, 10 ML MgO:N are 
deposited at Td = 300 K and pn 2 = 5 ' 10~ 6 mbar. No 
nitrogen was found in the sample, i.e. x < 0.2 %. Thus, 
if Mo is covered by MgO, the catalytic effect of the sub- 
strate is inhibited. 

Figure 4 compares STS curves of 11 ML films of 
MgO0.96N0.04 and undoped MgO. 18 Doping by N leads 
to a shift of the Fermi level Ep towards the valence band 
by about 1 eV and an additional peak with maximum 
at 0.3 eV below the conduction band minimum (CBM). 
This is observed for three different N concentrations, i.e. 
MgO thicknesses, as probed by STS. The peak energy 
varies laterally by ±0.3 eV. Density functional (DFT) 
calculations of bulk MgO with substitutional N predict 
occupied p-levels close to the valence band maximum 
(VBM) and an unoccupied, spin polarized level within 
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FIG. 4. (Color online) a) dI/dU(U) spectra measured by 
STS on a 11 ML high MgOo.9eNo.o4 island (U Bta , b = 3 V, 
/stab = 0.5 nA, [/ m od = 40 mV) at several positions (straight 
lines) and on a 11 ML thick pristine MgO island (dashed 
line); average film thickness in both cases: 7 ML. b) Calcu- 
lated charge density of the unoccupied N-induced states of 
a N-N dimer at the surface. Large green and small white 
spheres show Mg and O respectively, c) Calculated density 
of states (DOS) for N-N dimer at MgO surface, HOS: high- 
est occupied state; all states between -1.2 eV and 3 eV are 
N-induced; CBM is marked by arrow. 



the middle of the band gap, if self interaction correction 
is included. 11 ' 12 

This disagrees with our experiment. Since surfaces 
arc not included in these calculations, we performed 
first-principles DFT calculations including the surface 
within the spin-polarized generalized gradient approxi- 
mation 26 using projector augmented-wave potentials as 
implemented in the Vienna Ab initio Simulation Pack- 
age (VASP). 27 Correlation effects on the p-shclls of N 
dopants are accounted for by the DFT+U scheme in Du- 
darev's approach 28 with an on-site effective Coulomb pa- 
rameter U ff=3.4 eV. A kinetic energy cutoff of 500 eV 
and a 6 x 6 x 1 T-centered fc-point mesh was used. The 
supercell consists of 9 atomic layers MgO(OOl) using the 
experimental lattice parameter and a 16 A thick layer of 
vacuum. All atomic positions as well as the thickness of 
the MgO slab were fully relaxed. 

Three different configurations with N atoms in the sur- 
face layer have been calculated: (i) one N-atom substi- 
tuting an oxygen, (ii) one N-atom at the interstitial site 
and (iii) an N-N dimer with one N atom being substi- 
tutional and the other at the nearest interstitial site. In 
case (ii), the calculation was initiated with an N inter- 
stitial, but the system relaxed to a configuration where 
N and O exchanged their places, i.e. N ended up substi- 
tutionally and O intcrstitially. However, the N derived 



p-states for case (i) and (ii) are found within 2.2 eV 
above VBM very similar as for MgO bulk. 11 ' 12 We con- 
clude that a more complex N structure as, e.g., the N-N 
dimer at the surface shown in Fig. 4(b) and exhibiting 
unoccupied N-type p-states close to CBM (Fig. 4(c)), is 
responsible for the observed peak. 

In conclusion, we prepared thin films of N-doped MgO 
with an N-concentration up to 6 % by using the catalytic 
effect of Mo(OOl) for N2 dissociation. Compared with 
pristine MgO, an additional state close to the conduction 
band minimum has been observed by STS which could 
not be attributed to simple N impurity configurations. 
Helpful discussions with P. Mavropoulos and S. Parkin as 
well as financial support by SFB 917-A3 and HGF.YIG 
VH-NG-409 are gratefully acknowledged. 
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